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Abstract 

We present the hrst evidence for B semileptonic decays into the charmed baryon based on 
420 fb-i of data collected at the r(4S') resonance with the BABAR detector at the PEP-II e~^e 
storage rings. Events are tagged by fully reconstructing one of the B mesons in a hadronic decay 
mode. We measure the relative branching fraction 13{B Af X£~Vi)/B{B —>■ A'^ / A^ X) = (3.2 ± 
0.9stat. i 0.9syst.)%- The significance of the signal including the systematic uncertainty is 4.9 
standard deviations. 
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1 Introduction 



The B decays to charmed baryons are not as well understood as the decays into charmed mesons. In 
particular, there is limited knowledge, both theoretical and experimental, about the B semileptonic 
decays into charmed baryons. If the charmed baryonic production is dominated by the emission of 
an external W boson, as in the case of the mesonic production B DX, then we can expect the 
rate of the semileptonic events to be about the same for the baryonic and mesonic processes: 

B(B DXi-Vj) B(B A+Xe-Vj) 

B(B ^ D/DX) ^ B(B ^ At/A:X) ^ ' 

where £ = e or For the mesonic process, the ratio is B{B DXl^V()/B{B D/DX) k, 10% [1]. 
If the baryonic ratio is found to be smaller than the mesonic ratio, there is a significant contribution 
from internal W emission in the baryonic production. The Feynman diagrams of B decays with 
internal and external W emissions are shown in Fig. 1. 




Figure 1: B decays with an external W emission (left) and an internal W emission (right). 

About 90% of the inclusive semileptonic B X£~V£ branching fraction [2] can be accounted 
for by summing the branching fractions from exclusive B — >■ D^*\Tr)i~T'£ decays. Semileptonic 
B decays to charmed baryons may account for some of the remaining difference. These decays 
have not yet been observed. A previous search for B A^Xe~T'e by the CLEO collaboration [3] 
obtained an upper limit of B(B AfXe~Ue)/B(B A+ /A~ X) < 0.05 at the 90% confidence 
level. 

In this paper, we present the first evidence for B A'^Xe^Vc decays*^, where X can be any 
particle(s) from the B semileptonic decay other than the leptons and the A^. The B — > A'^Xe~T'e 
signal yield is obtained by a fit to the yl+ invariant mass. We perform a blind analysis using Monte- 
Carlo (MC) samples and the yl+ invariant mass sidebands on data to optimize the selection criteria 
and estimate the backgrounds. We also present the results for a similar search for B —>■ A'^ Xji'V^j^. 

2 The BABAR Detector and Dataset 

This analysis is based on data collected with the BABAR detector at the PEP-II storage rings. 
The total integrated luminosity of the dataset is 420 fb""^ collected on the T(4S') resonance. The 

^Charge-conjugate modes are implied throughout this paper. 
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corresponding number of produced BB pairs is roughly 460 million. An additional 42 fb~^ data 
sample taken at a center-of-mass (CM) energy 40 MeV below the T{4S) resonance is used to 
study background from e"*"e~ ff {f = u,d,s,c,T) events (continuum production). The BABAR 
detector is described in detail elsewhere [4]. Charged-particle trajectories are measured by a 5-layer 
double-sided silicon vertex tracker and a 40-layer drift chamber, both operating in a 1.5-T magnetic 
field. Charged-particle identification is provided by the average energy loss (dE/dx) in the tracking 
devices and by an internally reflecting ring-imaging Cherenkov detector. Photons are detected by 
a CsI(Tl) electromagnetic calorimeter. Muons are identified by the instrumented magnetic-flux 
return. A detailed GEANT4-based MC simulation [5] of BB and continuum events has been used 
to study the detector response, its acceptance, and to test the analysis techniques. 

3 Simulation of B ^ A'^Xi~Vi Decays 

Due to a lack of available theoretical models for semileptonic B decays to charmed baryons, we 
use an ad-hoc model to study selection optimization in our analysis. In our model, the B decays 
semileptonically into charmed baryons through an intermediate massive particle Y ss B ^ Y£~V£, 
according to a phase space model [6]. The Y subsequently decays into a A'^, an anti-nucleon 
(anti-proton or anti- neutron) , and ni (^2) charged (neutral) pions, with ni and n2 distributed as 
Poissonians with a common mean of 1.25. We apply isospin symmetry in the final states of the Y 
decay. 

Due to the ad-hoc nature of the signal model, we tune its parameters as part of the analysis. 
We perform a first optimization of the selection criteria using a pseudo-particle Y mass my of 4.0 
GeV/c^ and a width Ty of 0.6 GeV/c^, with the Y decay constrained by ni+n2 < 4. After unblind- 
ing the A'^ invariant mass signal region, we compare the sideband-subtracted signal distributions 
of the A'^ and lepton momentum spectra, and the charged and neutral pion multiplicity, with the 
corresponding ones from the signal MC simulation. We then tune the signal model parameters 
to resemble the observed distributions on data. For instance, by adjusting the mass and width of 
the pseudo-particle Y, we can control the shape of the lepton momentum spectrum. We find a 
good agreement with the electron spectrum in data for my = 4.5 GeV/cP, Ty = 0.2 GeV/c^, and 
n-i + n2 < 6. We use the rctuncd signal model to rc-optimizc the selection criteria and to esti- 
mate signal efficiency. We check a posteriori that the two-pass selection criteria optimization does 
not strongly depend on the initial signal model parameters used. The Af and lepton momentum 
spectra for a sample of signal MC events are shown in Fig. 2 for two the models before and after 
tuning. 

4 Event Selection 

We select B semileptonic decays in events containing a fully reconstructed B meson (Stag), which 
allows us to constrain the kinematics, to reduce the combinatorial background, and to determine 
the charge and flavor of the signal B, up to flavor-changing mixing. 

We first reconstruct the B semileptonic decay, selecting a lepton with momentum p| in the CM 
frame higher than 0.35 GeV/c. Electrons from photon conversion and tt^ Dalitz decay are removed 
using a dedicated algorithm, which uses a least-squares fit to reconstruct the vertex between two 
tracks of opposite charges whose kinematical parameters are compatible with a photon conversion 
or a TT^ Dalitz decay. Tracks identified as electrons or muons are required to be within 0.1 cm 
of the the interaction point at their point of closest approach in the transverse plane. Candidate 
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Figure 2: A'^ and lepton momentum spectra for a sample of signal MC events using the models 
before tuning (dashed line) and after tuning (solid line) as described in the text. 



baryons, with the correct charge-flavor correlation with the lepton, are reconstructed in the 
pK~Tr~^, pKg, pKg7:~^7r~ , A7r~^, Att^tt^tt^ modes. Kg (A) candidates are reconstructed in the 
TT+TT" {pTT~) decay mode. 

In events with multiple B — > A'^Xi^Vi candidates, the candidate with the highest A^-i~ vertex 
fit probability is selected. 

We reconstruct Stag decays of the type B DY', where Y' represents a collection of hadrons 
with a total charge of ±1, composed of n^vr^ + n2K^ + n'^Kg + u'^tt'^, where + 71-2 < 5, n'^ < 2, 
and 77,4 < 2. Using {D~^) and D*^ (D*^) as seeds for B^ (B^) decays, we reconstruct about 
1000 types of decay chains. For each of the i?tag decay modes, the purity V is estimated using MC 
simulation. V is defined as the ratio of signal over background events with uies ^ 5.27 GeV/c^. 

The kinematic consistency of a Stag candidate is checked using two variables: the beam-energy 

substituted mass ruES = y^/^ ~ Pb^ ™^ energy difference /S.E = Eb — \/s/2. Here ^/s refers 
to the total CM energy, while pB and Eb denote the momentum and energy of the -Stag candidate 
in the CM frame. For correctly identified -Btag decays, the tues distribution peaks at the B mass, 
while /^E peaks at zero. We select a -Btag candidate in the signal region defined as 5.27 GeV/c^ 
< rriEs < 5.29 GeV/c^, excluding Stag candidates with daughter particles in common with the 
charmed baryon or the lepton from the B semileptonic decay. 

In the case of multiple Stag candidates, we select the one with the largest purity V of the Stag 
mode; in the case of multiple candidates with the same Stag mode we select the one with the 
smallest \^E\ value. 

The Stag, -^t ^ s-^d ^~ candidates are required to have the correct charge-flavor correlation. 

By fully reconstructing one S in the event and by requiring the presence of a proton or A 
candidate, used for the A'^ reconstruction, the resulting sample shows a high purity, as is generally 
the case for a tagged analysis. In order to minimize model-dependent effects that can be introduced 
by a specific set of selection criteria, we only require the lepton momentum to be greater than 0.35 
GeV/c, as described above, and the missing momentum to be greater than 0.2 GeV/c. This last cut 
removes background from hadronic S ^t-^ decays in which all the particles in the X system 
have been reconstructed and one hadron is misidentified as a lepton. We also require the total 
charge of the reconstructed event to be zero, in order to reduce combinatorial background in the 
Stag reconstruction from missing particles. 

To obtain the S semileptonic signal yields, we perform a one-dimensional binned maximum 
likelihood fit to the A'^ invariant mass distribution. Backgrounds to the process of interest can 
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be divided according to whether they contain a correctly-reconstructed A'^ with a mass value in 
the signal region (peaking), and those that do not. MC studies of generic BB and continuum 
events show that the peaking background comes mainly from hadronic B decays, with a 

fake electron from gamma conversions or vr'^ Dalitz decays, or a hadron misidentified as a muon. 
The number of peaking background events from hadronic B A'^X decays is estimated from the 
simulation, as: 

iVpeak = 2 • epeak " B{B ^ yl+X) • iV^^ 

where e^eak is the efficiency of reconstructing fake B A'^Xi~V£ events in a hadronic B A'^X 
sample computed as the ratio of reconstructed and generated events, and is the number of 
BB pairs corresponding to the data luminosity. We use B{B A'^X) = (3.9 ± 1.2)%, which is 
the statistical average of the charged and neutral B branching fractions [1]. We estimate A'peak = 
5.2 lb l.Ostat. i 1.7syst. events and A'pcak = 15.3 ± 1.4stat. i 4.9syst. events for the electron and muon 
samples, respectively. The sources of systematic error are described in Sec. 6. 

The yl^ invariant mass distribution is fitted with the sum of three probability density functions 
(PDFs): one Gaussian for semileptonic B A'^Xi'Vg events, another Gaussian for peaking 
background events from hadronic B A'^X decays, and a first order polynomial for combinatorial 
BB and continuum background events. The two Gaussians share the same mean and width, and 
the width is fixed to the value obtained by fitting the hadronic B A'^X data sample, as described 
in Sec. 5. The amount of peaking background PDF is fixed to the MC prediction. We first fit the 
A'^ invariant mass sidebands, defined as the mass window of 2.23 — 2.26, 2.31 — 2.34 GeV/c^, to 
constrain the background PDF parameters. We then fit the A'^ invariant mass including the signal 
region, where the free parameters are the signal and background yields, and the A'^ mass mean 
value. 



5 Measurement of Branching Fractions 

The A'^ invariant mass distribution is compared with the results of the fit in Figure 3 for the 
B A'^Xa.'ui decays, where we show separately the electron and muon sample. 
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Figure 3: Fit to the A'^ distribution for B A^Xe^Ve (left) and B A^Xfi^V^ (right): the 
data (points with error bars) are compared to the results of the overall fit (solid line). The peaking 
background contribution is shown with a shaded area. The combinatorial BB and continuum 
background is shown with a dashed line. 

In Fig. 4 and 5, we show the distributions for the A'^ and electron momentum spectrum, 
and the charged and neutral pion multiplicity in the X system for the B A'^Xe~Ve sample. 
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These distributions are sideband-subtracted, by selecting events in the invariant mass signal 
region, and subtracting the combinatorial BB and continuum background using the invariant mass 
sidebands, and then corrected bin- by-bin by the efficiency estimated from the signal MC. The 
peaking background contribution is subtracted using the MC prediction. 
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Figure 4: Sideband-subtracted, efficiency-corrected, normalized distributions for the and elec- 
tron momentum spectrum m.B ^ A^Xe~V(, decays: the data (points with error bars) are compared 
to the signal MC prediction. 
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Figure 5: Sideband-subtracted, efficiency-corrected, normalized distributions for the charged and 
neutral pion multiplicity in the X system m. B ^ A^Xe~V(. decays: the data (points with error 
bars) are compared to the signal MC prediction. 

In order to reduce the systematic uncertainty, the exclusive B{B A'^ Xl^Vi) branching frac- 
tions are measured relative to the hadronic B{B — > A'^ / A^ X) branching fraction. To determine 
the hadronic branching fraction, we use a sample of B{B — > A'^ / A^ X) events selected similarly to 
the semileptonic channel. We select the reconstructed A'^ candidate with the highest vertex prob- 
ability, and a Stag candidate in the signal region defined as 5.27 GeV/c^ < niES < 5.29 GeV/c^, 
excluding Stag candidates with daughter particles in common with the charmed baryon. In the 
case of multiple Stag candidates, we select the one with the largest a priori purity of the Stag rnode; 
in the case of multiple candidates with the same Stag niode (same a priori purity), we select the 
one with the smallest \^E\ value. 

To obtain the B — > A'^ /^c ^ signal yield, we perform a one-dimensional binned maximum 
likelihood fit to the A'^ invariant mass distribution with the sum of two PDFs: a single Gaussian 
for hadronic B A^/A^X events, and a first order polynomial background for combinatorial 
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BB and continuum background events. The A'^ invariant mass distribution is compared with the 
results of the fit in Figure 6. 



^1400 
0I2OO 
glOOO 
d 800 
^ 600 
I 400 
W 200 



" 2.24 2.26 2.28 2.3 2.32 2.34 

A, Mass (GeV/c^) 

Figure 6: Fit to the A'^ distribution for B A'^/A^X: the data (points with error bars) are 
compared to the results of the overall fit (solid line). The combinatorial BB and continuum 
background is shown with a dashed line. 




Table 1: Signal yields and reconstruction efficiencies for the B A^Xi and B j A^ X 

decays with statistical errors. 

"Decay Mode TVdata e (xlQ-^) " 

B ^ A^Xe-V^ 38 ± 10 2.06 ±0.17 
B ^ AtXll-V^, 7.4 ±8.4 1.07 ±0.12 
1432 ±81 3.02 ±0.13 



The relative branching fraction B{B A'^Xi i'()/B{B j A^ X) is obtained by correcting 

the signal yields obtained from the fit by the reconstruction efficiency ratio: 

B{B ^ AtXl~V,) ^ /AUni\ / ehad \ 

B{B At/H^X) V ^^had / VesemiJ ' 

Here, A'^gemii (-^had) is the number oi B ^ ^t-^ signal events for the semileptonic (hadronic) mode, 
reported in Table 1 together with the corresponding reconstruction efficiencies e. 

6 Systematic Uncertainties 

By measuring the B — > A^X£~T'£ branching fraction relative to the B — > A'^ / A^ X branching 
fraction, many uncertainties in the reconstruction efficiency ratio cancel out. In particular, the 
uncertainties in the A'^ and -Stag reconstruction efficiencies and the vl^ decay branching fractions 
do not contribute. 

We categorize the remaining systematic uncertainties into additive uncertainties which directly 
affect the signal yield, and multiplicative uncertainties which affect only the branching fraction ratio. 
The systematic uncertainties that have been considered are described below and summarized in 
Tab. 2. 
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Table 2: Table of systematic uncertainties. The additive errors (errors on the signal yield) are 
shown in units of events. The multiplicative errors (errors on the reconstruction efficiency) are 
shown in units of percent. 



Additive systematics 


B A+Xe-V^ 


B ^ A+Xii-V^ 


Peaking background statistics (ev.) 


1.0 


1.4 


Uncertainty in B(B A'^X) (ev.) 


1.6 


4.7 


Lepton misidentification rate (ev.) 


0.7 


2.0 


Fit bias (ev.) 


0.3 


1.2 


Total additive; (c;v.) 


2.0 


5.1 


Multiplicative systematics 


B A+Xe-Ve 


B ^ A+Xfi-Vi, 


Reconstruction efficiency statistics (%) 


8.4 


11.4 


Model dependence (lepton momentum) (%) 


22.6 


71.8 


Model dependence {A'^ momentum) (%) 


7.5 


7.5 


Lepton identification efficiency (%) 


1.1 


2.7 


Selection order (%) 


3.7 


44.6 


Total multiplicative (%) 


25.5 


85.7 



6.1 Additive Systematics 

Systematic uncertainties in the signal yield are dominated by the peaking-background yield esti- 
mate. We evaluate this uncertainty by propagating the uncertainty in the B A'^X branching 
fraction, and the Poisson error due to the limited hadronic MC statistics. We also vary the lepton 
misidentification probabilities by 15% for both electrons and muons, and add in quadrature the 
corresponding variation in the peaking background rate. To evaluate the effect of a bias in the fit 
technique, we perform ensembles of MC experiments, in which events are generated according to 
the PDF shapes measured on data, varying the signal to background rate, and fitting for the signal 
as in the full analysis. The difference between the fitted value of the yield and the true value is 
taken as a systematic error. 

6.2 Multiplicative Systematics 

Systematic uncertainties which affect the reconstruction efficiency ratio are dominated by the uncer- 
tainty in the signal model. To estimate this, we look for deviations in the reconstruction efficiency 
as we vary the tuning parameters of the signal model. The deviation in the electron spectrum 
is taken as the Ax^/n.d.f. = 1 difference, where Ax^ is the data-MC difference divided by the 
statistical error on data, added in quadrature for each bin of momentum, and n.d.f. = 6 is the 
number of bins. The deviation in the A'^ spectrum is taken as the variation in the efficiency as 
a function of Af momentum. We also include the Poisson error contribution of the limited signal 
MC statistics. We estimate the systematic uncertainty due to particle identification by varying the 
electron (muon) identification efficiency by 2% (3%). Because the event selection order is slightly 
different in the semileptonic and hadronic sample selections, the reconstruction efficiency system- 
atics do not exactly cancel. We evaluate the corresponding systematic uncertainty by reversing the 
order of the lepton and Stag selection, taking the difference in the signal yield, corrected by the 
reconstruction efficiency, as the systematic error. The large systematic uncertainty coming from 
the selection order in the muon channel is due to the low statistics of the muon sample. 
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7 Results 



We measure the following branching ratios: 

B{B ^ a+/a;x) 

B{B ^ A+Xfj-Vf,) 
B{B ^ At/XX) 



(3.9±1.0stat.±l.lsyst.)% (2) 
(1.5±1.7stat. ±1.7syst.)%. (3) 



The result for the electron channel is compatible with the CLEO result [3] of ^ ^t-^^ ^e) 
/ B{B A'^/A^X) < 0.05. Despite the lower value of the muon result, the two channels are 
compatible within their errors. Taking the weighted average of the two channels, we obtain: 

^[^ ^ '^J'l = (3.2±0.W±0.Vt.)%. (4) 
BiB ^ At/A^ X) 

By using the hadronic branching fraction B{B A'^ / A^ X) = (4.5 ± 1.2)% [1], we obtain: 

B{B^AtXe-Ve) = (1.8±0.5stat. ±0.7syst.) X 10-=* (5) 
B{B^ A^Xirv^) = (6.6±7.5stat. ±7.6syst.) X 10-^ (6) 
B{B ^ Atxr-ui) = (1.5 ± 0.4stat. ± O.Gsyst.) x 10"^ (7) 

We evaluate the significance 5 of our result using the log likelihood ratio —2 log(>Cmax/'Co) where 
-Cmax is the maximum likelihood and Cq is the likelihood value fixing the signal yield to be zero, 
including statistical and systematic errors: 

Se = ^l-2\og{C^^/Co) = 4.6 (8) 



5^ = ^-21og(£„iax/>Co) = 0.8 (9) 

where the subscript e and /i denotes the electron and muon channels, respectively. The probability 
P that our measured signal results from a background fluctuation are computed from these S 
values to be Pe = 2.2 X 10~^ and = 0.22. In order to compute the significance for the combined 
result, we take the product of these two probabilities: P^ = 4.8 x 10~^. Converting this back to 
the significance, we obtain 5/ = 4.9 for the combined result. 

In conclusion, we present the first evidence for B semileptonic decays into the charmed baryon 
A'^ . The relative branching fraction B{B — A'^Xi~T'()IB{B — A'^ / A^ X) is found to be smaller 
than the corresponding one for the D charmed meson. We measure the A'^ and electron momentum 
spectrum in the B semileptonic decay, which could be of guidance to further development in the 
theoretical modeling of these decays. 
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